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The ground state properties of a one-dimensional system with particle-hole symmetry, consist- 
ing of a gate controlled dot coupled to an interacting reservoir, are explored using the numerical 
DMRG method. It was previously shown that the system's thermodynamic properties as a func- 
tion of the gate voltage in the Luttinger liquid phase are qualitatively similar to the behavior of a 
non-interacting wire with an effective (renormalized) dot-lead coupling. Here we examine the ther- 
modynamic properties of the wire in the charge density wave phase, and show that these properties 
behave quite differently. The number of electrons in the system remains constant as a function of 
the gate voltage, while the total energy becomes linear. Moreover, by tuning the gate voltage on the 
dot in the charge density wave phase it is possible to drive the wire through a first order quantum 
phase transition in which the population of each site in the wire is inverted. 
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The properties of one dimensional (ID) interacting 
systems have attracted much interest going back half a 
century 1,2 . Much recent effort has concentrated on un- 
derstanding the conductivity and I-V characteristics of 
a Luttinger liquid (LL) coupled to an impurity^. These 
properties are probed experimentally by measurements 
of the temperature dependent conduction through ID 
systems^, and tunneling spectroscopy into ID wires'. 
The measurements of Ref. |5| also indicate a localization 
transition in the wire for low densities which might be as- 
sociated with a charge density wave (CDW). Recently it 
has been realized that the occupation of an impurity level 
or a quantum dot coupled to an interacting ID sample 
may be used as a probe to its properties'^. The occu- 
pation is experimentally accessible, for example by the 
charging effect of the impurity on a quantum point con- 
tact in its vicinity 8 . A generic model for such a situation 
is a dot coupled to a lead. If such a dot is controlled by 
a gate one can change its orbital energy. In this paper 
we shall demonstrate that by following the occupation of 
the dot as a function of the gate voltage it is possible to 
identify a transition in the wire from a LL to a CDW 
phase. In the CDW regime of the wire the occupation of 
the dot switches abruptly once the dot orbital crosses the 
Fermi energy. Moreover, by changing the gate voltage of 
the quantum dot the wire is driven through a first order 
quantum phase transition (QPT), in which the ground 
state of the entire system switches. 

The difference between measuring the conductivity 
through the dot-lead system or the local density of states 
at the impurity 9 , and probing the dot occupation using 
e.g. a quantum point contact must be emphasized. It 
is commonly known that the special behavior of the LL 
spectrum near the Fermi energy causes the conductivity 
to vanish, even in the presence of a single impurity 10 . 
Thus measuring a conductivity of a system containing 
a LL lead and a dot is not expected to show any level 
broadening at zero temperature. 

Nevertheless, connecting the dot to the LL lead, while 
measuring its occupation using a quantum point contact, 



results in a different situation. In this case the LL acts 
as a reservoir for the dot, and a level broadening can 
be seen. This difference results from the fact that such 
a level is coupled to all the lead states, even those far 
from the Fermi energy, and thus the special behavior of 
the LL near Ep is not expected to dominate. It was in- 
deed shown, that the dot's level width follows the regular 
Breit-Wigner line shape, even in the presence of interac- 
tions, although the strength of coupling between the dot 
and the ID wire is renormalized as a function of the in- 
teractions in the wire£. 

In this paper we show that a qualitatively different be- 
havior of the dot's population is exhibited once the wire 
is in the CDW phase. It is known that a one-dimensional 
system of spinless electrons with nearest neighbor repul- 
sive interactions undergoes a QPT between a LL phase 
for weak interactions, and a CDW phase for strong 
interaction o 11 ^ 2 ' 13 . The different behavior of the dot 
population between the LL and the CDW phases stems 
from the difference between the ground states of the two 
phases in the wire. It is argued that while the LL is only 
locally influenced by the coupling to the dot, a semi- 
infinite CDW coupled to the dot undergoes a QPT when 
the orbital's level crosses the wire's chemical potential. 
We show that a sharp jump in the dot's population oc- 
curring at this point is accompanied by an abrupt switch 
of the ground state of the entire system. This switch is 
not a simple level crossing, but rather a true QPT since 
it becomes a sharp transition only once the wire is semi- 
infinite 14 . Other hallmarks of a first order QPT, such as 
a discontinuity in the first derivative of the grand poten- 
tial, as well as an inversion of the CDW order parameter 
as a result of the inversion of the wire's occupation for 
each site, are seen. 

The Hamiltonian describing a spinless fcrmionic sys- 
tem composed of a single orbital ("dot") coupled to a 
ID nearest neighbor interacting wire ("lead"), is given 
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by H = H d ot + Hdot-iead + Hlead, where 

Hdot = e a*a, (1) 
Hdot-iead = -V(a)c\ + c\a), 

L-l 

Hlead = -tyi(c]Cj+l +h.C.) 
3=1 

L-l x l 

+ i^i^Cj- -){c] +1 c ]+1 - -). 

3=1 

We denote the dot's energy level by eo, V (t) is the dot- 
lead (lead) hopping matrix element, and I is the nearest- 
neighbor interaction strength in the lead, a' (a) is the 
creation (annihilation) operator of an electron in the dot, 
and Cj (cj) is the creation (annihilation) operator of an 
electron at site j in the lead. In the interaction term, a 
positive background is assumed. The lead hopping ma- 
trix element, t, is taken as 1, in order to set the energy 
scale. 

In order not to increase the number of free parameters, 
the model discussed here does not contain a dot-lead in- 
teraction term. Nevertheless, we have checked that such 
a term does not influence qualitatively the results of this 
paper. We refer the reader to a previous work£, in which 
it was explicitly shown that an addition of such an in- 
teraction term influences only the effective position and 
width of the dot level. 

As was previously mentioned, while the lead is in the 
LL phase, the Breit-Wigner formula with an effective dot- 
lead coupling describes the dot population n(eo) quite 
welt. For a dot coupled to a LL lead in which the inter- 
action parameter g < 1/2, a jump in the dot population 
was predicted by Furusaki and Matveev 6 . However, for 
the spinless model with nearest neighbor interaction the 
LL parameter is always in the range of 1/2 < g < 1, and 
the Furusaki-Matveev jump is not expected. 

The grand potential Cl = H — /j,N (N is the parti- 
cle number operator) was diagonalized using a finite-size 
DMRG calculation^^ for different values of V and /, 
with a lead of up to L = 300 sites. In order to obtain 
particle-hole symmetry, for which the LL to CDW tran- 
sition was widely discusse d 11 ' 12 ! 13 , \i is set to zero. The 
total number of particles in the system is not fixed during 
the renormalization process, so that the results describe 
the experimental situation of a finite section of a ID wire 
which is coupled to a dot and to an external electron 
reservoir. In particular, the total occupation obtained 
for such a system can be non-integral. As a function of 
the dot's level (eo), the following ground state properties 
were calculated: the system's grand potential f2, the total 
number of electrons N and the dot population ridot ■ The 
population of lead sites were also calculated in order to 
differentiate between local effects of the dot population 
and global phenomena in the lead. 

We begin by comparing the behavior of n ( iot in the two 
phases of the interacting lead. For spinless fermions with 
nearest-neighbor interactions, it is known that the system 



experiences a phase transition (of a Kosterlitz-Thouless 
typei&), between the LL and the CDW phases, at I = 2t. 
Indeed, the results presented in Fig. Q] show a qualita- 
tive difference between ndoti^o) in these two regimes. In 
the LL phase the curves fit quite well the Breit-Wigner 
formula- with an effective coupling V e f / 

«dot(^,eo)= (2) 



T J-2t ^-(l--^r) + ((l-^-)e-e )2 

while in the CDW phase the width becomes zero, and a 
jump in ridot occurs at e = 0. This jump is associated 
with the degeneracy of the CDW ground state. 
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FIG. 1: The dot population ridot as a function of the level 
energy eo- For I < 2 (LL phase) the curves fit the non- 
interacting formula with an effective coupling constant V e // 
(lines - DMRG results, symbols - fit to Eq.[2j), while for I > 2 
(CDW phase) the width is zero. Exactly at the transition 
point (/ = 2) Udot is continuous but doesn't fit the non- 
interacting formula. 

In order to support this argument, we first calculate 
the properties of the ground state for a lead not coupled 
to a dot, by taking V = t and the limit eo - * 0. The 
CDW order parameter can be defined asii 

P(i) = i(-iy[2n« - n(i - 1) - n(i + 1)], (3) 

where n(i) =< eje^ > is the occupation of the i'th lead 
site in the ground state. The value of P(i) does not 
change much as a function of the spatial coordinate, ex- 
cept in the vicinity of the lead's edges. We thus define 
P = P{L/2), 

In Fig. [2] the dependence of P on the interaction 
strength /, as obtained by the DMRG method, is shown, 
and compared to the exact results^. For I < 2, the sys- 
tem is in the LL phase and indeed P = ± 10" 4 . In this 
case the population of each lead site is 1/2 and the lead 
is half filled; i.e. N = L/2. 

As expected, the metal-insulator transition occurs at 
1 = 2. For / > 2, in which the system is in the CDW 
phase, the value of P is finite. For values of I which are 
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FIG. 2: The CDW order parameter P as a function of the 
interaction I for a 300-sites lead (DMRG - symbols, theory - 
lines). For the CDW phase P is inverted between the cases 
eo — * + and eo — * _ , while in the LL phase both cases 
result in P — 0. Inset: P as a function of eo for I = 1.5 
(triangles), 2.5 (squares) and 3.5 (circles). The LL case results 
in a constant P = 0, while an inversion of P occurs at eo = 
for the two CDW cases. 



far from the transition point (7 > 2.5), we get a very good 
fit to the theory. Near the transition point the numerics 
tend to emphasize the charge oscillations, resulting in 
too large values of \P\. This tendency, however, does not 
affect the following qualitative conclusions. 

As can be seen, the CDW ground state is different for 
the cases eo — * + and e — * 0~, resulting in two values 
of P (i.e., ±|P|) for any value of I in the CDW regime. 
For eo = the ground state is two-fold degenerate in the 
thermodynamic limit. 

In the CDW case, special care should be devoted to the 
number of electrons in the system, and to the difference 
between even and odd lead length L. For eo — > and 
V = t, the lead length is effectively L + 1. Denoting the 
ground state for eo — > ± by ipQ, an d it's population in 
each site as nf z , where j = for the dot and 1 < j < L 
for the lead sites, the particle hole symmetry implies that 
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= 1 — everywhere. For odd L (even L + 1) one 



has the requirement rij ' = n^lj, an d the number of 

electrons in both i/>o +) and ip^ is N = (L + l)/2. When 
L is even (L + 1 is odd), however, one has an additional 

symmetry requirement, n 
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tpQ ' , and the states differ in their electrons number: ip, 
contains N = L/2 electrons, for any I > 2, while ipQ 
has N = L/2 + 1. It is worth noting that n{ +) > 1/2, 
while < 1/2. 

When eo ^ 0, i.e., a dot with a finite on-site energy 
is connected to the lead, it has some local influence on 
the ground state in its vicinity. For a lead in the CDW 
phase, the two states described above are slightly modi- 
fied, and can be denoted by ip^ and tp^~K Nevertheless, 
the main influence of the dot is to lift the degeneracy be- 
tween those states. If eo < 0, the dot state population 
is high and is energetically preferable due to the 
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dot-lead hopping term. For eo > the opposite happens, 
resulting in a preference of The occupation of the 

dot-lead system in the CDW phase is shown in Fig. [3] 
When the dot's orbital energy changes from a negative 
to a positive value, the system switches from to 
ip( + \ In the following we show, by a calculation of the 
free energy, that this is indeed a first order QPT. The 
resulting phase diagram is drawn in the inset of Fig. 3] 




FIG. 3: The occupation of the dot and the first lead sites of 



the CDW states tpg +> (circles) and (squares), of a lead 

with 300 sites and 1 — 4 which is coupled to the dot (filled 
symbols) with V = 0.8. For eo = these two states are 
degenerate in the thermodynamic limit. Lines are guide to 
the eye. 



Before discussing the free energy in the CDW phase, 
we turn to an exact calculation of the energy and of the 
total number of electrons in a non-interacting case. Using 
Green functions technique^, one can find that the change 
in the total density of states in the system, due to the 
presence of the dot, is 



A^(e) = is|-ln(e-eo-£(e)), 
7T oe 



(4) 



where S(e) = (V/t) 2 e/2 + i(V 2 /t)y/l - (e/2t) 2 is the self 
energy of an electron in the do1j£. 

Therefore, the change in the number of electrons in the 
entire system is 



AiV( M ,e ) 



Au(e)de 



(5) 



2/ 



11 _ x (j, - e - SRS(/x) 
2 + ^ tan 3S(^ ■ 



and the change in the free energy of the system is 
AO(/i, eo) = AE — nANQi) = (6) 
{e-n)Av{e)de = - AN(e,e )de. 

-2t J-2t 

In Fig. [4] the total number of electrons in the system 
is presented. For I = 0, N fits the predicted formula 
(Eq. [5]) well. As noted above, the agreement between 
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the numerical results for a finite lead and the exact re- 
sult for a semi-infinite lead, is obtained due to the fact 
that the number of particles can vary during the DMRG 
process. For the LL phase (0 < I < 2), N(eo) looks 
quite similar, varying between L/2 + 1 (at eo — > — oo) 
and L/2 (at eo — -> oo), taking the average value (L + l)/2 
at eo = 0. As in the case of ndot, the results fit Eq. © 
with a renormalized dot-lead coupling for small values of 
I. Increasing I towards the transition point (7 > 1.5), 
results in a less accurate fit. 



ible. Hence, infinitesimal changes of the electrons num- 
ber are possible. On the other hand, the CDW phase 
is insulating and thus incompressible, which results in a 
constant N. 
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FIG. 4: The total number of electrons in the system, as a 
function of eo, for L = 299, y, = 0, V = 0.3, and different 
values of /. Symbols - DMRG results. Lines - fit to Eq. [5] 
Inset: the phase diagram: a Kosterlitz Thouless transition 
(solid line) occurs at I = 2, and a first order phase transition 
(dashed line) at eo = for I > 2. 

The values of V e ff(I) obtained by the fit of the N 
curves to Eq. ((SJ), are in good agreement with the values 
obtained by the fit of ndot to Eq. fl^} . We find that V e j / 
decreases monotonically with increasing /, exhibiting the 
RPA like behavior described in Ref. [Z|. 

The CDW phase (I > 2), however, is qualitatively 
different: for an odd lead length A(eo) remains a constant 
integer (N = (L + l)/2) which does not depend on eo at 
all. For even L, N = L/2 + 1 for e < 0, and N = 
L/2 for eo > 0. This is a direct result from the switch 
of the system ground state from tp(> to at eo = 

0. For an odd lead length, the total number of sites 
in the system (L + 1) is even, so that N is equal for 
both states. For an even lead length, L + 1 is odd, and 
the total number of electrons is changed by one when 
eo passes 0. Except for the decrease of one electron at 
eo = for an even lead, N remains constant in the CDW 
phase, independent of eo. Thus even with the continuous 
change in the population of the dot as a function of the 
dot's level for eo ^ 0, the number of electrons in the entire 
system remains constant. The change in the occupation 
of the dot as a function of eo is compensated by the lead. 

This difference in the behavior of A as a function of 
eo between the two phases (i.e., constant for the CDW 
phase, compared to a continuous decrease for the LL 
phase) is a direct manifestation of their transport prop- 
erties. The LL phase is metallic, and therefore compress- 




FIG. 5: The ground state free energy as a function of eo 
for different interaction strengths. The lead-dot coupling was 
taken as V = 0.3, and the interaction strength I takes values 
between and 3.5, in jumps of 0.5. The lines represent the 
DMRG results (dashed line -1 = 2), while the symbols show 
the fit to Eq. [6] for values of / lower than 2. Inset: fi(eo) — 
fi(— eo) as a function of eo for / = 0,1,2,3 (symbols) and 
fi(eo) — fl(— eo) = eo (line). For all values of interaction we 
get |fi(eo)-«(-eo)-eo| < 10~ 6 . 

In Fig. [5l typical numerical results for the free energy 
O(eo) — O(0) as a function of eo are shown. A perfect fit 
between the DMRG results and the exact formula Eq. © 
for / = is obtained. In the LL phase (0 < I < 2) our 
DMRG calculations show that the effect of interactions 
on fl can be fitted by replacing V in Eq. ^ by the same 
effective coupling V e f / obtained for the behavior of ndot 
discussed in Ref. 0, and for the behavior of N discussed 
above. For the CDW phase (I > 2), however, there is 
obviously a qualitative change in the energy curve: the 
dependence of Q on eo is linear both below and above 
e = 0, with an abrupt change of ^ at e = 0. 

These results point out that the single impurity, con- 
nected at one end of a long interacting lead, has a well de- 
fined influence on the ground state of the entire coupled 
system. As discussed above, when the dot level passes 
through eo = 0, the lead's population is inverted at every 
site, leading to an inversion of the CDW order parame- 
ter, as presented in Fig. [2^inset). The dot population is 
inverted as well. As a result, a dramatic change in the 
dependence of the free energy of the system on eo occurs. 

Since is the free energy of the system, the jump in 
it's first derivative might be a sign of a first order QPT. 
In order to see whether this non-analyticity of the free 
energy is just a trivial level crossing (LC) of two levels 
in the system, the dependence of the transition shape 
on the system size L is explored. For the case that the 
sharp transition in the energy results from the fact that 
the external field (in our case the gate voltage) commutes 
with the Hamiltonian, and thus a LC is possible, no size 
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FIG. 6: The ground state free energy as a function of eo 
for different CDW system sizes, with V = 0.3 and / = 2.5 
(dotted, dashed-dotted and dashed lines). Inset: the results 
for L = 29 are compared to the case of / = —2.5, which is 
inside the ferromagnetic phase, and to the CDW case with 
L = 30. 



dependence of the sharpness of the transition is expected. 
On the other hand, a real QPT will become sharp only 
in the thermodynamic limit (i.e., semi infinite lead). 

In order to compare the two scenarios (LC vs. QPT) 
we solve the same Hamiltonian (Eq. (I)) with strong 
attractive interactions, i.e. I < —It. In this case, 
transforming the Hamiltonian to an XXZ spin chain, the 
strong attractive interactions regime corresponds to the 
ferromagnetic phased (while the CDW is equivalent to 
the Nccl phase). In the fermionic case these attractive in- 
teractions yield a two-fold degenerate ground state, com- 
posed of states in which the sites are either entirely occu- 
pied or empty (in case there is an explicit restriction to 
half filling, a "phase separation" evolves). Pinning by the 
dot lifts the degeneracy, because an occupied dot causes 
a preference of an empty lead, and vice versa. Thus the 
influence of the dot on the coupled wire is superficially 
similar to the CDW case in so far as in both cases the 
ground state of the entire system is determined by the 
dot orbital. Nevertheless, the ferromagnetic system is 
clearly a LC system, and no dependence on the system 
size is expected. 

The dependence of f2 on eo was calculated for L = 
29, 99, 299 in the CDW case, and for L = 29 in the ferro- 
magnetic case. From the results shown in Fig. [5] (and its 
inset) it is clear that these two cases are different. The 
ferromagnetic system is indeed a trivial LC system, with 
two competing states whose energies cross each other for 
eo = 0. Since these two states are eigenstates of the 
Hamiltonian even for a finite L, the size does not play a 
role, so that even for L — 29 one can see a sharp transi- 
tion between the two ground states. 

On the other hand, for a small (but with an odd L) 
CDW system f2 shows a smooth dependence on eo j with- 
out any non-analyticity. As a matter of fact, for any finite 



system the electron levels are expected to be mixed, re- 
sulting in avoided crossing of the two lowest many-body 
levels^. In other words, the CDW states are not true 
eigenstates of the Hamiltonian for a finite system. In- 
deed, indications for the sudden jump of ^p- are seen 
only for larger system sizes (i.e., for L > 200). 

The dot-dashed line in Fig. [SJinset) shows that for a 
short lead of an even size (L = 30 in that case) there is, 
however, a non-analyticity in f2 as a function of eo. The 
comparison between even and odd lead sizes emphasizes 
our conclusion stated above. While for an odd lead size 
ip(~' and tfj^ have the same number of electrons, for an 
even size of the lead, the transition from ?/>( _ ) to i/>( + ' at 
eo = involves a decrease of the electrons number by one. 
The Hamiltonian H conserves the number of particles, 
so that in the case of L even, these two states are not 
coupled, and the transition between them is a simple LC, 
not showing a size dependence. For an odd L, however, 
tp^ and ip^ are coupled by H, so that for a finite 
L they are actually mixtures of the CDW states, thus 
presenting sharper r2( e o) dependence for larger systems. 
We thus conclude, based on these two comparisons, that 
in the case of a CDW with an odd L this transition is a 
QPT, which happens for L — > oo, i.e., for a semi-infinite 
learJ . 




FIG. 7: The dependence of the dot occupation on eo for differ- 
ent CDW system sizes, with V = 0.3 and I — 2.5. Inset: the 
absolute value of the slope in the limit eo — > as a function 
of L for short leads, together with a linear fit. 

As pointed out above, in the numerical results for large 
lead sizes, the non-analyticity of the free energy is seen 
clearly, so that it is obviously a first order transition. 
However, a similar conclusion can be drawn from the 
results of short leads, by scaling the results against the 
lead size L. Although for finite size systems the order 
parameter changes continuously, the slope of this change 
grows with the system size L. For a first order transition 
in d dimensions one expects a power law dependence of 
the slope as L d , while for a second order transition the 
power law should be fixed by some universal exponents^. 
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For the CDW model, the population of the dot plays 
the role of an order parameter, and one can check its 
behavior near eo = for different lead sizes. In Fig.[7]the 
occupation as a function of eo is shown for some short 
lead sizes (L = 29,49,69) and for a large one (L — 299). 
It is indeed seen that while for a long enough lead there 
is a jump in the dot population, for short leads there is a 
continuous change in the occupation, with a larger slope 
for larger lengths. Zooming into the regime of eo = we 
find that although L > 200 is required in order to see 
a clear discontinuity in il, L > 100 is long enough for 
showing a jump in ndot- The dependence of the slope on 
L, for L < 100, is shown in the inset of Fig. [7] As one 
expects from finite size scaling predictions for first order 
phase transition, there is a very good linear fit of the the 
order parameter on the lead size. 

We now turn to another consequence of the particle- 
hole symmetry. For eo ^ the symmetry results in 
the fact that the ground states for +eo and for — eo 
(■0( + ) and respectively) have an inverted popula- 

tion everywhere. It is thus clear that (ip( + ' |-ff; ea( i|V^ + ' ) ) = 
(¥- ) \Hlead\¥- ) ), and similarly (<#<■+) \H dot -iead\^ +) ) = 
(^iHdot-ieadl^)- This implies 



Q(e ) - fi(-eo) 
(^\H dot \^) 



(7) 



e on^ ot - (-eoKTot = £ o, 



where nj ot 



( n dot) represents the dot populations of i/>( + ) 
The last equality results from the symmetric pop- 



ulation n dot = 1 



b dof 



The relation Eq.[7]does not depend on the interaction, 
and should exist for both the LL and the CDW phases. 
Obviously it should be obeyed for the non-interacting 
case, where the energy is given by Eq. ([6|), thus leading 
to the following nontrivial relation 



tan de 



¥1 
i 



eo, 



(8) 



which is evident since the LHS is the occupation of the 
dot£ when [i > 2t. 

In Fig. [5jinset) a plot of f2(eo) — fi(— eo) as a function 
of eo is shown, for values of / between and 3. As can 
be seen, Eq. ([7]) is valid for all values of I. 

In conclusion, we have shown that the occupation of 
a dot coupled to a one dimensional wire can be used to 
identify the different phases in the wire. It is known 
that as a function of the interaction strength / the wire 
goes through three different phases (which may also be 
mapped onto the phases of an XXZ spin chain). For 
/ < — 2t the ground state is doubly degenerate, with all 
sites equally full or empty, corresponding to the ferro- 
magnetic state of a spin chain. In the intermediate range 
of interactions (— 2i < I < 2t) the translational symme- 
try is unbroken and a Luttinger Liquid phase occurs. For 
I > 2t the ground state is doubly degenerate correspond- 
ing to the two possible CDW states or the two different 
antiferromagnetic realizations for the spin chain. 

For the Luttinger liquid phase the occupation of the 
impurity does not show a jump when the impurity level 
crosses the Fermi energy. In the two other phases the 
impurity level splits the degeneracy by favoring one of 
the two ground states depending on whether the impu- 
rity level is empty or filled. Nevertheless, we have shown 
that the physics in these cases is different. While for 
the ferromagnetic phase a simple level crossing occurs 
with a sharp jump in occupation of the impurity for any 
length of the wire, in the CDW phase the position of 
the impurity level drives a first order QPT in the ther- 
modynamical limit between the two CDW states, while 
for a finite wire the jump is smeared. This phase tran- 
sition shows all the hallmarks of a first order QPT, such 
as a size dependence, a jump in the order parameter and 
a discontinuity of the derivative of the grand canonical 
potential. 
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in which the result of the integral on the LHS does not 
depend on the parameters V and t. 

A physical insight into Eq. ||SJ) can be gained by taking 
the derivative of both sides with respect to eo, and using 
the definition of the self energy given above. This yields 



r-21 



3E(e) 



W-« (e-e -KS(e))2 + (3<S(e)): 



:de = 1, 



(9) 
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